Hawking radiation as seen by an infalling observer 
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We investigate an important question of Hawking-like radiation as seen by an infalling observer 
during gravitational collapse. Using the functional Schrodinger formalism we are able to probe the 
time dependent regime which is out of the reach of the standard approximations like the Bogolyubov 
method. We calculate the occupation number of particles whose frequencies are measured in the 
proper time of an infalling observer in two crucially different space-time foliations: Schwarzschild 
and Eddington-Finkelstein. We demonstrate that the distribution in Schwarzschild reference frame 
is not quite thermal, though it becomes thermal once the horizon is crossed. We approximately fit 
the temperature and find that the local temperature increases as the horizon is approached, and 
diverges exactly at the horizon. In Eddington-Finkelstein reference frame the temperature at the 
horizon is finite, since the observer in that frame is not accelerated. These results are in agreement 
with what is generically expected in the absence of backreaction. We also discuss some subtleties 
related to the physical interpretation of the infinite local temperature in Schwarzschild reference 
frame. 



I. INTRODUCTION 

Understanding quantum radiation from collapsing ob- 
jects, despite tremendous progress made in the last thirty 
years, still represents a major challenge in theoretical 
physics. An asymptotic observer, watching gravitational 
collapse of a massive object, will start registering radi- 
ation of quanta coming from the fields excited by the 
non-trivial metric of the background space-time. As the 
radiation progresses, it acquires more and more thermal 
features. Finally, when the horizon is formed, radiation 
becomes completely thermal, in agreement with the fact 
that radiation from a pre-existing horizon must be ther- 
mal, as shown for the first time by Hawking [lj. In this 
picture, radiation is treated in the context of quantum 
mechanics, while the background is described by classi- 
cal general relativity. Moreover, the background is held 
fixed, i.e. the form of the metric and even the mass of 
the collapsing object (or a black hole) are held fixed. It 
is extremely difficult to include any form of backreaction 
in the realistic models (see @, H, 0, [E Hj 0] for some 
attempts) . 

In the absence of any backreaction, the radiation flux 
as seen by an asymptotic observer is constant at late 
times. If integrated over the infinite time (an asymptotic 
observer sees formation of a horizon only after infinite 
amount of his time) this flux diverges. However, in clas- 
sical general relativity, an infalling observer (say an ob- 
server which is falling together with a collapsing object) 
will cross the Schwarzschild radius in the finite amount of 
his proper time. This time interval corresponds to an in- 
finite interval of the asymptotic observer. The question 
then is what such an observer would see. The answer 
will likely strongly depend on the space-time foliation 
one chooses. If an observer is infalling in Schwarzschild 
reference frame, his trajectory in any given moment cor- 
responds to an accelerated observer. Though such an 
observer reaches the Schwarzschild radius in finite time, 
his acceleration diverges there and the temperature he 
measures will diverge. Such an observer would need to 



encounter all of the radiation that an asymptotic observer 
would see, only in finite time. Will such an observer be 
burned by an intensive radiation as he is sailing toward 
the Schwarzschild radius? Not necessarily. A "particle" 
has no invariant meaning. The distribution of produced 
particles (occupation number as a function of the fre- 
quency) strongly depends on the frequency of produced 
" particles" . Since proper time coordinates for an asymp- 
totic and infalling observers are different, it is likely that 
these two different observers will register different distri- 
butions of particles. 

A very important complementary picture is given by 
Eddington-Finkelstein coordinates where the singularity 
at the horizon is absent. Such an observer would still 
register radiation due to the time dependent metric, how- 
ever we do not expect that such radiation will be ther- 
mal in the whole frequency range. In addition, since the 
Schwarzschild radius is not a singular point, we do not 
expect the temperature of such radiation to diverge there. 

Unfortunately, arguments of this type have remained 
only qualitative so far. The reason is that the notion of 
particles is well defined only in the asymptotically flat 
regions (e.g Minkowski or Schwarzschild). In the Bo- 
golyubov method, which is widely used in this context, 
one matches the coefficients between the two asymptot- 
ically flat spaces, i.e. Minkowski at the beginning and 
Schwarzschild at the end of the gravitational collapse. 
The mismatch of these two vacua gives the number of 
produced particles. What happens in between is beyond 
the scope of the Bogolyubov method. Thus, the question 
what an infalling observer (in different space-time folia- 
tions) would register during the collapse can not be an- 
swered in the context of the Bogolyubov method. How- 
ever, recently developed functional Schrodinger formal- 
ism, goes beyond the approximations of the Bogolyubov 
method. In this context, one can successfully attack the 
question of radiation as seen by an infalling observer, 
which is the goal of this paper. 
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II. RADIATION AS SEEN BY AN INFALLING 
OBSERVER IN SCHWARZSCHILD REFERENCE 
FRAME 

We consider a spherical domain wall representing a 
spherical shell of collapsing matter. The wall is described 
by only the radial degree of freedom, R(t). The metric is 
taken to be the solution of Einstein equations for a spher- 
ical domain wall. The metric is Schwarzschild outside the 
wall, as follows from spherical symmetry 8| 



d S 2 = -(l-^^+a-^r 1 ^ 2 

r r 



2 dtt 2 , r > R(t) 
(1) 

where R s = 2GM is the Schwarzschild radius in terms of 
the mass, M, of the wall, and 



d9 2 + sin 2 



(2) 



In the interior of the spherical domain wall, the line ele- 
ment is flat, as expected by Birkhoff's theorem, 



ds z 



-dT z + dr z 



r 2 d6 2 + r 2 sin 2 Odcj) 2 



r < R(t) (3) 



The equation of the wall is r = R(t). The interior time 
coordinate, T, is related to the asymptotic observer time 
coordinate, t, via the proper time of an observer moving 
with the shell, r. The relations are 



dT 



1/2 



and 



dt 



l 

B 



B 



dR 



1/2 



where 



R 



(4) 



(5) 



(G) 



We will single out an observer who is falling together 
with the collapsing shell, which means that he moves on 
the trajectory r — R(t) and measures the proper time r. 
By integrating the equations of motion for the spherical 
domain wall, Ipser and Sikivie 0] found that the mass is 
a constant of motion and is given by 

1 , 



M = -[Vl + i? 2 + yjB + i? 2 ]47TCTi? 2 (7) 

where R T = dR/dr, while a is the surface tension (en- 
ergy density per unit area) of the wall. It is assumed 
that max(i?) < (47rGcr) _1 to avoid the case in which the 
domain wall is already within its own Schwarzschild ra- 
dius to begin with. This expression for M is implicit 
since R s = 2GM occurs in B. Solving for M explicitly 
in terms of R T gives 



M = AiraR 2 ^! + R 2 T - 2-kGuR]. 



(8) 



From the mass, the velocity is given by 



\Rr\ 



\ AiraR 2 



2-koGR I 



1. 



(9) 



Now consider a massless scalar field $ which propa- 
gates in the background of the collapsing shell. The ac- 
tion for the scalar field is 



S 



(10) 



where g^ v is the background metric given by Eq. (TTJ) and 
Eq. ([3]). Decomposing the (spherically symmetric) scalar 
field into a complete set of real basis functions denoted 
by {fk(r)} 



$ = ^a fc (r)/ fe (r) 



[11) 



we can find a complete set of independent eigenmodes 
{bk} for which the Hamiltonian is a sum of terms. The 
total wavefunction then factorizes and can be found by 
solving a time-dependent Schrodinger equation of just 
one variable. 

Since the metric inside and outside the shell have dif- 
ferent forms, we can split the action (fTU|) into two parts 



,fl(r) 

dr J drr 2 



Vl + R 2 



{d T <S>) 2 



-y/rn%{d r *)' 



(12) 



So 



2tt / dr 



drr 2 



R{r) 



B 



y/B + B*l-Rs/r 



j£±3 (i -£)(».)■ 



(13) 



where we used ((!]) and ((5]). 

The most interesting things happen when the shell 
(and the infalling observer sitting on the shell) ap- 
proaches the Schwarzschild radius. From Eq. ([9]) we see 
that R T is constant in the limit when R — > R s . There- 
fore the kinetic term for Si n is roughly constant. The 
kinetic term in S ou t goes to zero as R — > R s , so the Si„ 
kinetic term is dominant. Similarly the potential term 
in Si n goes to a constant while the potential term in 
Sout becomes very large, so the potential term in S ou t 
dominates. Strictly speaking, this argument fails in the 
neighborhood of r = R, however the dominant contribu- 
tion to the integrals will be for r ^ R. Therefore we can 
write the action as 



S« 2tt 



dr 



1 



drr 2 — 

Rs 



'drr 2 ^ f 1 



B 



(14) 
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where we have changed the limits of integration from 
R(t) to R s since this is the region of interest. 

Using the expansion in the modes Eq. (fTTj) . we can 
rewrite the action as 



S 



dr 

\_Rr_ 

2B 



1 



1 



ak(T)C kkl a k '(T 



dfc(r)A fefe /d fe '(T) 



(15) 



where a = da/dr, and A and C are matrices that are 
independent of R(t) and are given by 



A fcfc < =4tt / drr z f k {r)f k ,{r) 
h 



C kk , = 4tt 



(16) 



(17) 



From the action (|5ip we can find the Hamiltonian, and 
according to the standard quantization procedure, the 
wave function ip(a k , r) must satisfy 



(18) 



dtp 



(19) 



2 

where 



i^i + ^n^A-^fcfc-nit 



\Rr\ 

2B 



a k {j)C kk ia k >{T) 



n. 



da k {r) 



(20) 



is the momentum operator conjugate to a k (r). 

So the problem of radiation from the collapsing domain 
wall for the infalling observer is equivalent to the problem 
of an infinite set of coupled harmonic oscillators with 
time dependent frequency. Since A and C are symmetric 
and real (i.e. Hermitian), it is possible to simultaneously 
diagonalize them using the principal axis transformation. 
Then for a single eigenmode, the Schrodinger equation 
takes the form 



1 rrrm d , l*M 



Kb 2 



i>(b,T) 



.dip(b, 



ch 



(21) 

where m and K denote eigenvalues of A and C, and b is 
the eigenmode. 

Re-writing Eq. (|2ip in the standard form we obtain 



*(M)=<^ (22) 



1 d 2 m o, , l2 



where 



K 



\Rr\ 



dr) 
\R T \ 



B^l + R 2 By/1 + R% 



(23) 



and 



'/ 



J dr'^/T+Rj. 



(24) 



where we defined w2 = K/m. The exact solution to 
Eq. (H]) is given by |g| 



7T/CT 



1/4 



exp 



TL ( 2jl + — U 2 

2 \ p p 2 ' 



. (25) 

where p^ = dpj dr\ and p is given by the real solution of 
the ordinary (though non-linear) differential equation 



p m + cj 2 (?y)p = i 



with initial conditions 



p(0) = — , p„(0) = 0. 



The phase a is given by 
a(r?) = 



(26) 



(27) 



(28) 



Complete information about the radiation in the back- 
ground of the collapsing shell is contained in the wave- 
function l[2"5"|). 

Consider an observer with detectors that are designed 
to register particles of different frequencies for the free 
field $. Such an observer will interpret the wavefunction 
of a given mode b at some later time in terms of simple 
harmonic oscillator states, {if n }, at the final frequency, 
to. The initial (r = 0) vacuum state for the modes is the 
simple harmonic oscillator ground state 



^) = ^V /4 e -™^/2 



(29) 



The relationship between the coordinates 77 and r is given 
by Eq. (|2~4"|) . which near R s simplifies to drj/dr =const. 
The number of quanta in eigenmode b can be evaluated 
by decomposing Eq. (|2"5|) in terms of the states and eval- 
uating the occupation number of that mode. The wave- 
function for a given mode in terms of simple harmonic 
oscillator basis is given by 



where 



db<p* n (b)^{b,v) 



(30) 



(31) 



which is an overlap of the wavefunction at some later 
time tp(b, r) with the simple harmonic oscillator basis 
functions. The occupation number at eigenfrequency ui 
is given by 



(32) 
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The occupation number in the cigenmode b is then given 
by (see Appendix A) 



N(t,lo) 



top 

where p T = dp/ dr. 

N 



\ 2 / \ 2 
1 \ (p. 



up' 



ioR s = 2,A, 10 



vjp 



(33) 
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FIG. 1: The occupation number iV as a function of the proper 
time t/Rs for various fixed values of particle frequencies ujR a . 
The curves are lower for higher values of ujR b . The occupa- 
tion number diverges as the infalling observer approaches R s , 
which happens as r — » r c = 2. 33333. ..-R^ 1 ■ 

For fixed u>, N is a function of time since p and 
p T are functions of time. In Fig. Q] we plot the occu- 
pation number of produced particles as a function of 
time (for several fixed frequencies QR S ). The amount 
of proper time needed for the shell (and the infalling 
observer) to reach R s can be obtained by integrating 
Eq. ([9]). For a = 0.01-Rj 3 this critical proper time is 
t c = 2. 33333. --Rj 1 . Fig.[T]shows that, as the infalling ob- 
server approaches R s , the occupation number increases 
and diverges exactly at R s . The same conclusion can be 
drawn by analyzing the occupation number N in Eq. (|3"3")) 
as a function of p and p T (see Appendix B). This is in 
agreement with what one would expect in the absence 
of backreaction. Hawking showed in [l[ that, for an 
asymptotic observer, the flux of particles at late times 
is steady (constant in time). This means that it diverges 
for a fixed background (i.e. fixed mass of the collaps- 
ing object) since we effectively keep pumping energy into 
the collapsing object so that its mass is constant despite 
radiation. For the asymptotic observer it takes infinite 
amount of his time for the collapsing object to reach its 
own Schwarzschild radius. This infinite interval corre- 
sponds to a finite time interval of the infalling observer's 
time. Thus, one may conclude that the infalling observer 
has to encounter an infinite number of particles before 
he hits R s [l(|. However, we should keep in mind that 
it is the particle occupation number that is divergent, 
not the actual number of particles detected by the in- 
falling observer. The number of particles detected may 



be quite different from the occupation number. As the 
infalling observer approaches R s , he and his detector get 
blue-shifted. His detector can not register particles whose 
wavelength is larger than his detector; for them the de- 
tector is "inside" the particles [lfj. Thus, an infalling 
observer might not be able to register all of the created 
particles. We only calculated the occupation number of 
particles contained in the wave function describing radi- 
ation at some later time. The occupation number does 
not tell us where the particle are nor how they propagate 
further. 

t/R s = 2.33, 2.333, 2.3333 




FIG. 2: The occupation number N as a function of fre- 
quency ujRs for fixed values of the proper time r/R 3 . The 
occupation number increases for larger values of t/R s as 
r -> r c = 2.33333.. .R' 1 . 

For fixed time r, the occupation number N is a func- 
tion of frequency Q at some fixed r. From Eq. (|23[) we 
see that, in order to keep the values of Co fixed in time, 
luq — > as B — > 0. Thus, Co varies with too and not with 
time since all of the values for Co must be calculated at the 
same final time. In Fig. [5] we plot the occupation number 
N as a function of frequency CoR s for fixed values of the 
proper time t/R s . We can compare these plots with the 
occupation numbers for the thermal Planck distribution 



Planck 



i 



1 



(34) 



where j3 is the inverse temperature. The curves have 
manifest non-thermal features — the occupation number 
does not diverge at CoR s = as in the Planck distri- 
bution and small scale oscillations which were absent in 
the Planck distribution are present here. However, as 
t — > t c = 2. 33333. --Rj 1 the distribution becomes more 
thermal. 

In Fig. [3] we fit the best fit linear approximation of 
the spectrum in Fig. [51 From Eq. I|34p we see that the 
slope of the line is f3, which is proportional to the inverse 
temperature, i.e. T ex Several important features 

of the Hawking- like radiation can be read off this plot. 
Non-thermal features of the radiation are apparent. It 
is interesting that, exactly as argued in (Tlj . departures 
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t/R s = 2.333, 2.3333, 2.33333 

Ln(l + 1/N) 




calculated is the stress-energy tensor near the horizon. 
As pointed out in [l3T ] , where a simple 1 + 1 dimensional 
model was studied, the local value of the stress-energy 
tensor is small since the local vacuum polarization can- 
cels out the divergent temperature energy density due 
to radiation. However, one should keep in mind that a 
complete realistic 3 + 1 dimensional calculation does not 
exist. Also, it may happen that a local backreaction is 
small but a global integrated backreaction can not be ne- 
glected since a collapsing object excites fields which are 
not localized in space and time. These are the questions 
of great importance in black hole physics but are outside 
of the scope of this work. 



FIG. 3: Plot of ln(l + 1/N) as a function of frequency ujR 3 
for fixed values of the proper time t/R 3 . The slope of the 
best fit line is (3, which is the inverse temperature. The 
non-thermal features disappear and the temperature diverges 
as the Schwarzschild radius is approached, i.e. r — > r c = 
2.33333.. .R7 1 . 



from thcrmality are larger for larger frequencies. As the 
time progresses and the infalling observer is approach- 
ing i? s , radiation becomes more and more thermal even 
at large frequencies. Finally, at R s it becomes purely 
thermal. At that point a black hole is formed and it 
is natural that the spectrum become thermal, as known 
from various studies of quantum radiation from a pre- 
existing horizon. We can not extrapolate the results for 
any further evolution of the infalling observer since the 
near-horizon approximation used in Eq. (|5ip fails. In 
the standard picture, the infalling observer will reach the 
singularity, i.e. a region of infinitely strong gravitational 
field, in some finite proper time, and any results about 
produced particles in that regime are likely to be non- 
physical. This may not happen if the quantum effects 
are capable of erasing the classical singularity at the cen- 
ter, just as in atomic physics quantum mechanics got rid 
of the singularity of the Coulomb potential which has an 
identical 1/r behavior (l2l |. 

It is apparent that the slope of the ln(l + 1/N) versus 
u)R s curve is decreasing as the infalling observer is ap- 
proaching R s . Exactly at R s , the slope of the curve is 
zero indicating that the temperature of the radiation is 
infinite. This is not surprising since, as it is well known, 
the asymptotic observer in the nearly flat asymptotic re- 
gion will register Hawking radiation with some finite tem- 
perature. When this temperature is blue-shifted back to 
i? s , it clearly diverges. Thus, the occupation number of 
particles as seen by an infalling observer will have a dis- 
tribution with the temperature which diverges as R s is 
approached. 

Infinite near horizon local temperature of the Hawking- 
like radiation may or may not indicate that the backre- 
action due to radiation will be important in that region. 
The standard lore is that such backreaction is small. For 
that question, the relevant quantity which needs to be 



III. RADIATION IN INGOING 
EDDINGTON-FINKELSTEIN COORDINATES 

Now we consider the collapse from the point of view 
of an ingoing Eddington-Finkelstein observer. This is a 
different space-time foliation than that in Schwarzschild 
coordinates, and we expect crucially different results. In 
particular, since the metric is not divergent at the hori- 
zon, we do not expect infinite temperature there. 

For this purpose, we define the ingoing null coordinate 
v as 

v = t + r* (35) 

where r* is the tortoise coordinate. We can then rewrite 
Eq. (I) as 



ds A 



R, 



dv 2 + 2dvdr + r 2 dQ 2 , r > R(v). 



(36) 

where the trajectory of the collapsing wall is r = R(v). 
The interior metric is the same as in Eq. (TSJ) The interior 
time coordinate, T, is related to the ingoing null coordi- 
nate, v, via the proper time on the shell, r. [Note that 
the proper time r is different from the same quantity in 
Schwarzschild coordinates since the space-time foliation 
is different.] The relations are 



dT / fdRY 

W 1+ UJ 



and 



where 



dr B \ dr V V dr 



R 



(37) 



(38) 



(39) 



Consider again a massless scalar field <I> which prop- 
agates in the background of the collapsing shell. The 
action for the scalar field is 



5* 



(40) 
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where g^ v is the background metric given by Eqs. ([3]) and 
([55)) . Decomposing the (spherically symmetric) scalar 
field into a complete set of real basis functions denoted 
by {fk(r)} 



(41) 



we can find a complete set of independent eigenmodes 
{b k } for which the Hamiltonian is a sum of terms. 

Since the metric inside and outside of the shell have 
different forms, we again split the action Eq. (|40|) into 
two parts 



S in = 27T i dT 



R(v) 



drr 2 [-(d T $) 2 + (d,$) 2 ] , (42) 



S out = 2ir I dv I drr" 



R, 



l - — (d r <f) 2 . (43) 



We are again interested in the near horizon behavior of 
the radiation, i.e. as R — > R s . In this limit we can write 
Eq. ([38]) as 



dv 



1 

2R~r 



(44) 



where R T = dR/dr. The explicit functional dependence 
of R T can be found again from Eq. which is just 
the consequence of the fact that M in Eq. ([7]) is the 
conserved quantity also in Eddington-Finkclstcin coor- 
dinates. Then with the help of Eq. (|37p we can write 
Eq. (g2]) as 



Si', 



rR(v) 

2ir I dv / drr 2 



1 



1 



2 ^R v /2(R v /2 + 1) V 

+ 2^R v /2(R v /2 + l)(d r $) 2 

(45) 



where _R„ = dR/dv. Obviously, the action is not singu- 
lar as R(v) — > R s , unlike the Schwarzschild case. From 
Eqs. ([33]) and (05]) we can write the total action as 



S «2tt / dv 



1 



1 



drr 2 — 

2 ^/R v /2(R v /2 +1) 



drr^^^O,.* + / drr 2 d r $d v ® 

JRs 

R s 



drr z 1 



(d r $)< 



(46) 



where we have changed the limits of integration from 
R(v) to R s since this is the region of interest. 



Now using the expansion in modes Eq. ([4"Tj) . we can 
rewrite the action as 



S 



dv 



1 



1 



:d k A kk > d k > 



2 y /R v /2(R v /2 + l 
-d k Y kk >a k / + -a k Y kk \a k i + -a k C kk ia k > 



(47) 



where d = da/dv, and A, Y and C are matrices that are 
independent of R(v) and are given by 



A kk , = 2tt 
Y kk , = Air 



C kk , = 8ir J 



R, 

r/ drr 2 f k (r)f k ,(r), (48) 
o 



drr 2 f k (r)f' k/ (r), (49) 



Rs 



f' k (r)f' k ,(r). (50) 



However if we take that the matrices are symmetric and 
real, we can see that Y = Y -1 , so we can write the action 
as 



S: 



dv 



1 



1 



d k A kk >d k ' 



2^/2(^/2 + 1) 

+ ^Yfcfc' («fe fl fe' + a kd k >) + -a k C kk ,a kl 

"(51) 

From the action Eq. (j5ip we can find the Hamiltonian, 
and according to the standard quantization procedure, 
the wave function ip(a k ,v) must satisfy 



W IT, 

ov 



(52) 



or 

.dip r 1 

dv 



^R v /2{R v /2 + ljn^A-^kfc/n* 



a k 



VR v /2(R v /2 + l)Y 2 kk ,(A )kk' + C kk i \ a k 



VR v /2(R v /2+l)U k Y kk ,(A- 1 ) kk , 



a,,. 



where 



da k 



(53) 



(54) 



is the momentum operator conjugate to a k . Using the 
momentum , we can rewrite the Schrodinger equation 
as 

.dip \l 
dv 



■^Rj2lRj2 + V ) U k {A- 1 ) kk ,ll k 



a k ^R, v /2(R v /2 + l)Y 2 kk , (A" 1 )**/ + C kk , a k 



ip 



i-^/Rj2jR~j2+T)Y kk ,{A~ 1 ) kk ,S kk , 



(55) 
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where Skk' is the Kroneckcr delta function. 

So the problem of radiation from the collapsing do- 
main wall for the infalling observer is equivalent to the 
problem of an infinite set of decoupled damped harmonic 
oscillators with time-dependent frequency. Since A, Y 
and C are symmetric and real, it is possible to simulta- 
neously diagonalizc them using the principle axis trans- 
formation. Then for a single eigenmode, the Schrodingcr 
equation takes the form 



dv 



d 2 



y/R v /2(R v /2 + 1) 



K b 



2m 



(56) 



where m, y and K denote eigenvalues of A, Y and C, 
and b is the eigenmode. 

Re-writing Eq. (|56|) in the standard form we obtain 



2m db 2 



2m, 



ij)(b, rj) 



i« (57 ) 
or) 



where 



fa) = —2 



K 



m 2 m ^R v /2{R v /2 + 1) 



vfi 2 yjR v /2{R v /2 + l) 



and 



7] = / dv'^/R v /2(R v /2 + l) 



(58) 



(59) 



where we defined u>q = K/m. To find solutions to equa- 



(60) 



tion Eq. ([57]) we use the ansatz 

VKM) = e- OT/2 ">(M). 

This leads to the equation for </>(&, rj) 



1 d 2 4> 
2m db 2 



moo 



ocp = i^—. 



2 " r drj 

As discussed in Ref. [|| , this has the implicit solution 

1/4 



(61) 



(p(b, rj) 



o ia (v) 



7Tp 



cxp 



vm / 
2 [p 



(62) 

where p^ is the derivative of the function p(rf) with re- 
spect to m and the defining equation for p is 



PTln +u 2 (r/)p 



1 



(63) 



The initial conditions for p are taken at some large value 
of rj (i.e. large value of R) denoted by rji, so that 



1 



, PriiVi) = 0. 



(64) 



The phase a is defined by 



1 f v 
«(»/) = — 2 J 



drj 
P 2 (jl') 



Then Eq. dBDJ give 



-yri/2m 



(65) 



(66) 



where ^> given in Eq. (|62p . 

Consider an observer with detectors that are designed 
to register particles of different frequencies for the free 
field $ at early times. Such an observer will interpret the 
wavefunction of a given mode b at late times in terms of 
simple harmonic oscillator states, {v? n }, at the final fre- 
quency to. The number of quanta in eigenmode b can be 
evaluated by decomposing the wavefunction Eq. (|66|) in 
terms of the states, {<p ra }, and by evaluating the occupa- 
tion number of that mode. To implement this evaluation, 
we start by writing the wavefunction for a given mode at 
time v < Vf in terms of the simple harmonic oscillator 
basis at v = vq 



where 



ip(b,v) = ^c n (v)ip n (b) 



dbtp* n (b)ij(b,v) 



(67) 



(68) 



which is the Gaussian overlap with the simple harmonic 
oscillator basis functions. The occupation number at 
cigcnfrcquency Q by the time v < Vf, is given by the 
expectation value 



(69) 



We evaluate the occupation number in the eigenmode b 
to be 



N{v,u>) = ^e^'/' 



1 - 



1 

Cup 2 



Cop 



(70) 



for v < Vf. 

In Fig. [4] we plot iV versus v/R s for various fixed val- 
ues of CoR s . The frequency Co is measured in time v. 
We can see that the occupation number at any frequency 
increases as v/R s decreases. Thus more particles are cre- 
ated as the shell reaches and crosses the horizon. How- 
ever, the number of created particles does not diverge as 
R(v) -> R s - 

We then numerically evaluate the spectrum of mode 
occupation numbers at any finite time and show the re- 
sults in Fig. [5] for several values of v/R s . The first sign 
of non-thermality is the fact that the occupation num- 
ber is non-divergent at Co = 0, as opposed to the thermal 
Planck distribution in Eq. (f3"4")l . 

For the values of parameters taken for plots in Fig. 0] 
and Fig. the Schwarzschild radius is crossed at v = 0, 
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coR s = 0.15, 0.2, 0.5, 1 

N 




FIG. 4: Here we plot JV versus v/R s for various fixed values of 
uiRs- The occupation number at all frequencies grows as the 
collapse progresses (i.e. v/R 3 decreases) but never diverges. 
The occupation number would probably diverge when the sin- 
gularity is hit at v = —0.126 because of the infinitely strong 
gravitational field, but our approximation breaks far from the 
Schwarzschild radius which is crossed at v = 0, so we do not 
extend our plots all the way to v — —0.126. The curves are 
lower for higher values of ljR s . 



v/R s = -0.01, 0.05, 0.1 

N 




FIG. 5: Here we plot N versus u)R a for various fixed values of 
v/Rs- The occupation number at any frequency grows as the 
collapse progresses (i.e. v/R s decreases) but never diverges. 

while the singularity is reached at v = —0.126. However, 
we can not extend our plots significantly beyond v = 
since there our approximation breaks. 

In Fig. [6] we plot ln(l + l/N) versus CuR s for various 
values oiv/R s . As v/R s decreases (as the shell is collaps- 
ing), the curves decrease. A thermal spectrum should 
gives us a straight line, however, we see that is not the 
case here. The best one can do is to fit the low fre- 
quency part of the spectrum and get the temperature 
in that regime. We see that the temperature (of the 
low frequency part of the spectrum) practically remains 
constant near the Schwarzschild radius (near v = 0). 
The numerical value that we get for the temperature is 



T w 0.7 /R s . Unlike the case of Schwarzschild coordi- 
nates, where the spectrum becomes thermal in the whole 
frequency range at late times, in Eddington-Finkelstein 
coordinates the spectrum never becomes thermal in the 
high frequency range. We can not extend our quanti- 
tative analysis all the way to the singularity, since the 
approximation we used in Eq. (|44[) fails far from R ss R s . 

We also note that all of the plots were made for the 
numerical value of the eigenvalue y defined after Eq. (|5rj|) 
of y — 1. Numerical experiments indicate that the spec- 
trum and the temperature do not change significantly for 
different values of the eigenvalue y. 

v/R s = -0.01, 0.001, 0.01, 0.05 

Ln(l + 1/N) 




FIG. 6: Here we plot \n(l + l/N) versus u>R a for various fixed 
values of v/R a . Clearly, one can not fit straight lines through 
these curves in the whole frequency regime. However, the 
low frequency part (w < Rg~ ) does appear thermal with the 
temperature that is practically constant in time v. We get 
the numerical value of T ~ 0.7/ R 3 



IV. CONCLUSIONS 

We studied the question of quantum radiation from 
a collapsing object from a point of view of an infalling 
observer. Precisely, we calculated the occupation num- 
ber of particles whose frequencies arc measured in the 
proper time of an infalling observer. As the collapsing 
object approaches its own Schwarzschild radius it excites 
fields that are propagating in this dynamical background. 
In the context of the functional Schrodingcr formalism, 
the information about the excited radiation is contained 
in the time dependent wave function. The distribution 
function of the radiation depends on the wave function 
overlap between the initial state (vacuum) and the state 
at some later time. This formalism allows us to track 
the time evolution of the radiation distribution function. 
This was not possible in approximations usually used in 
similar setups, e.g. in Bogolyubov method where one 
matches the coefficients between the two static asymp- 
totically fiat spaces, i.e. Minkowski at the beginning and 
Schwarzschild at the end of the gravitational collapse. 

We demonstrated several interesting results. Quantum 
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radiation accompanies gravitational collapse since the 
metric describing the collapse is time dependent. How- 
ever, the characteristics of the emitted radiation strongly 
depend on the space-time foliation of an observer. For the 
infalling observer in Schwarzschild coordinates the radi- 
ation distribution function is not quite thermal, though 
it becomes thermal when the collapsing object reaches 
its own Schwarzschild radius. We call such radiation 
Hawking-like or pre-Hawking radiation (as opposed to 
thermal Hawking radiation from a pre-existing horizon). 
In the absence of any backreaction, an asymptotic ob- 
server will observe a divergent flux of particles at infin- 
ity. However, the number of particles and their ener- 
gies have different meaning for different observers. The 
radiation distribution function (i.e. the particle occu- 
pation number as a function of frequency) depends on 
the time coordinate that the observer is using. Naively, 
an infalling observer would need to encounter an infinite 
number of particles, but in finite amount of his proper 
time. However, it is only the occupation number of parti- 
cles that is divergent, not the actual number of detected 
particles. An infalling observer can not detect particles 
whose wavelength is larger than his detector and there- 
fore would miss most of them. 



By fitting the approximate temperature corresponding 
to the radiation distribution function, we showed that 
the local temperature near the Schwarzschild radius, R s , 
is large and diverges exactly at R s . This is expected 
since this local temperature has to be infinitely resiftcd 
in order to get a finite temperature of radiation as seen 
by the asymptotic observer. After all, an infalling ob- 
server in Schwarzschild coordinates is accelerated at any 
given moment, with acceleration which diverge exactly 
at R s . Such an accelerated observer should "see" diver- 
gent temperature at the horizon in agreement with the 
Unruh effect [l4| . This still does not mean that the hori- 
zon is a violent place for the infalling observer since local 
vacuum polarization energy density is usually thought 
to cancel out the divergence of the temperature energy 
density, which yields a finite total energy density. 



APPENDIX A: NUMBER OF PARTICLES 
PRODUCED AS A FUNCTION OF TIME 

We use the simple harmonic oscillator basis states but 
at a frequency ui to keep track of different w's in the cal- 
culation. To evaluate the occupation numbers at time 
t < rj, we need only to set uj = (j-^t/), where the sub- 
script "f ' stands for "final" . So 



<p(b) 



fmu>\ 



1/4 e -mu>b 2 /2 



(Al) 



where H n are the Hermite polynomials. Then Eq. ([25 
together with Eq. (|31~j) gives 



1 



uir 2 p 2 



uir 2 p 2 



1/4 



V2 n n\ 



dCe- p < 2 / 2 H n (() 



1/4 



\Jl n n\ 



In 



where 



P= 1- 



v \ P P 2 



(A2) 



(A3) 



To find /„ consider the corresponding integral over the 
generating function for the Hermite polynomials 




- Z 2 (l-2/P) 



Since 



-z 2 +2z( 



n=0 



d" 



z=0 



(A4) 



(A5) 



(A6) 



Therefore 



The situation in Eddington-Finkelstcin coordinates is 
quite different. This coordinate system is not singular at 
the horizon, and an infalling observer in these coordinates 
is not accelerated. We find that the distribution of the 
particle occupation number is not thermal in the whole 
frequency range. By fitting the temperature only in the 
low frequency range (ui < Rj 1 ) we find the temperature 
T w 0.7/ Rs- This temperature is somewhat greater than 
the Hawking temperature as measured at infinity by an 
asymptotic observer = 1/(4ttR s ), which is in agree- 
ment with results obtained in different approaches (see 
e.g. [IH). However, the spectrum never becomes ther- 
mal in the whole frequency range, so it is difficult to talk 
about the temperature as seen by such an observer. 




Since 



H n (0) = (- 



(-l)"/ 2 V2^! (n ~l )!! 



(A7) 



(A8) 



and H„ (0) = for odd n, we find the coefficient c„ for 
even values of n, 




{Qp 2 ) 1 ^ V P 
For odd n, c„ = 0. 
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Next we find the number of particles produced. Let 



after rescaling can be written as 



2 

1 -P 



Then 



N(r,Q)= J2 



fQf?\P\ d£ 
2 



y 



n— even 

d l 



2 

" ^\p\(i~er /2 ' 

Inserting the expressions for £ and P, leads to 



1 - 



Cop 



(All) 



(A12) 



APPENDIX B: BEHAVIOR OF THE FUNCTION 
p NEAR THE SCHWARZSCHILD RADIUS 

To get an understanding of the number of particles cre- 
ated in the near horizon limit we need to investigate the 
behavior of the function p in Schwarzschild coordinates 
near the Schwarzschild radius. 

Near the horizon we can then write the velocity term 

as 

\Rr\ w const = A. (Bl) 

In this limit the position of the shell is then, from Eq. ([9]), 

R(t) w Ro - At (B2) 

where Rq is the initial position of the shell, so we can 
write 



Therefore the frequency becomes 



CB' 



(B3) 



(B4) 



Therefore the auxiliary equation becomes 



P ^ + UJ °C((Ro-R s )-ATy- P 3 
or using Eq. (|24[) we can write this as, 



Rs 



1 d 2 p 2 

C^d^ +UJ °C{{Ro-R s )-AT) r ^p 



(A10) fL = - 

dr' 2 



1 -T> 



1 

p 



(B5) 



where r' = Ar/(R - R s ), and / = ^(R S /C(R Q - 
R s )) 4 / 4 p. The boundary conditions are then 



/(0) 



Rs 



C(Rq — R s ) 



1/4 



df(0) 
dr, 1 



(B6) 



The last term with the l// 3 becomes singular as / — > 0. 
Let us consider another equation with this term replaced 
by something more well behaved. For example consider, 



dr' 2 



a 2 4 ,2 rV a c^ 



(R - Rs) 11 /* 
with boundary conditions 



.9(0) = 



Rs 



C(Rq — R s ) 



1/4 



dg(0) 

dr' 



0. 



(B7) 



(B8) 



Eq. (|B7|) implies that g(r') is monotonically decreasing 
as long as g(r') > 0. Furthermore, it is decreasing faster 
than the solution for / as long as / < 1, since the l// 3 
term in Eq. (|B5[) is a larger "repulsive" force than the g 
term in Eq. (|B7|) . Therefore we have 



Z(r') > ff(r') 



(B9) 



for all t' such that g(r') > 0. 

The solution for g is positive for all r 1 and, in particu- 
lar, 0(1) > for all the values of A 2 cj 3/2 C 5 / 4 i?s /4 /(i?o - 
R s ) n / 4 that we checked. Therefore f(r') is positive, at 
least for a wide range. 

Let fi = /(l) ^ 0. Then the equation for / can be 
expanded near t' = 1. 



d 2 / 
dr' 2 



(i2o - i?,) 11 / 4 



h 



1-t' 



1 

7? 



Integrating Eq. (|BT0|) we can then write 



dr/' 



,4 2 u; 3/2 ^ /4 C 5 / 4 



(iZo - Rs) 11 ' 



TUT-fi ln(l - r 



(B10) 



(Bll) 



as t' ^ 1. Hence /j(t = (i?o — R s )/A) is strictly positive 
and finite while p r (r = (i?o ~~ Rs)/A = — oo for finite and 
non-zero wo- 

We are calculating the occupation number N as a func- 
tion of frequency to at some fixed time. From Eq. (|B4p 
we see that, in order to keep u> fixed in time, coq — > as 
_B — * 0. Thus, w varies with wo and not with time. Since 
/ = (R S /C(R Q - R s )) 1/4 , and / -> (R S /C(R - R s )) 1/4 
for ojo — > 0, we see that p — > oo and p T — > as ojo — > 0. 
This implies that the occupation number N in Eq. 
diverges asnr c since B — > as r — > r c . 
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